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A transverse, comparatively low-frequency standing wave has been identified on the plasma column 
of a magnetically confined argon arc. The frequency of this wave, under some conditions, obeys a dis- 
persion relation given by Weinstock for a generalized hydromagnetic wave in an inhomogeneous, ro- 
tating plasma. The wave frequency, which may be higher than the ion gyrofrequency, is closely 
coupled with the fluid rotation velocity of the plasma. When the differences between these passes 
the ion gyrofrequency, transition to another, as yet unidentified, mode takes place. 



1. Introduction 

Abundant experimental observations on hydromag- 
netic waves in homogeneous, quiescent laboratory 
plasmas exist [see, for example, Jephcott, 1959; Nagao 
and Sato, 1960; Wilcox, DeSilva, and Cooper, 1961; 
Jephcott and Malein, 1964]. It is the purpose of this 
paper to present experimental observations on a low- 
frequency hydromagnetic wave in an inhomogeneous, 
rotating plasma. Such waves are qualitatively dif- 
ferent from the well-known homogeneous plasma 
waves in several important respects. First, the fre- 
quency of this wave is closely coupled, but not equal, 
to the frequency of the fluid rotation of the plasma. 
Second, the frequency of this wave, even when 
larger than the ion gyrofrequency, may be much less 
than that of an Alfven wave for the same k in a homo- 
geneous plasma (if an Alfven wave could exist under 
these conditions). Third, the frequency is much less 
dependent on the density of the plasma than on the 
magnetic field. The wave will be identified by com- 
paring its properties with those predicted by a disper- 
sion equation derived for a theoretical model of the 
rotating cylindrical plasma by Weinstock [1964]. 

The plasma used in this investigation is that of the 
magnetically confined, low pressure argon arc 
[McLane, Takeda, Thomas, and Thompson, 1965]. 
For the operation of this arc, three electrodes alined 
with a uniform (± 3 percent) magnetic field are utilized. 
Between the first, which is a 1/8 in. diam. thoriated 
tungsten cathode, and the second, which is a tungsten 
nozzle with a 1/8 in. bore, a 1/2 in., 80 A arc operates 
in argon at about 1 torr pressure. This arc acts as 
a virtual cathode for the main discharge, which 



operates at 8 A (for this work) in argon at 10 -3 torr 
between the nozzle and the third electrode, which is 
a hollow anode consisting of a stack of spaced tantalum 
washers. The length of the main discharge is variable, 
is referred to as the arc length, and denoted by L in 
what follows. Further experimental details on the 
discharge and on the method of measuring electron 
density and temperature distributions may be found 
in the reference cited above. 2 

The plasma of the main discharge is a very stable, 
cylindrical ribbon whose peak (axial) density is from 
10 19 to 10 20 m -3 and whose half density radius is from 
0.4 to 0.5 cm. The axial temperature is from 1.5 
to 3 eV and the half temperature radius is from 0.5 
to 0.9 cm. At the densities and temperatures quoted, 
the collision frequencies for momentum exchange may 
be compared as follows: the electron-ion collision 
frequency (10 8 — 10 9 sec -1 ) much exceeds the electron- 
neutral collision frequency (10 3 — 10 4 sec -1 ); the ion- 
neutral collisions (10 5 sec -1 ) dominate the ion-electron 
collisions (10 3 — 10 4 sec -1 ). The degree of ionization 
is from about 25 percent to over 75 percent. The 
ion temperature has been estimated by the same 
method used on a similar plasma by Lidsky et al. 
[1962] to be approximately 1 eV (error may be a factor 
of two either way). 

What concerns us here is that this plasma, under 
nearly all conditions of its existence, is in continuous 
oscillation. [McLane, Thomas, and Thompson, 1963]. 
The most commonly encountered modes are trans- 
verse, and it is our purpose to demonstrate that one of 
these modes is, in all probability, associated with a 
generalized hydromagnetic wave. 



1 On leave from the Nagoya University, Nagoya, Japan. 



2 It should be noted that the density and temperature distribution in the plasma, and the 
oscillation modes which are excited, will be influenced by the bore size and bore length of 
the nozzle electrode, but not appreciably by the current in the arc operating between 
cathode and nozzle. 



389 



2. Experimental Measurement of Oscillation 
and Rotation Frequencies 

Because it does not disturb the plasma, the best 
method of observing the oscillations is to use a lens to 
image the plasma column on a thin metal screen and 
place a photomultiplier behind a small hole in the 
screen at the place where the oscillation is to be ob- 
served. Similar results can be obtained with, in most 
cases, remarkably little disturbance to the plasma, by 
placing a small Langmuir probe (1/4 mm diameter X 1/2 
to 1 mm long) near it. The oscillations can then be 
observed in either the ion saturation current, or, better, 
in the floating potential. The photoelectric method 
was used for most of this work. 

Transverse waves of odd azimuthal wave number 
(m in the wave phase factor e i{(t)t+kz+me) ) are easily dis- 
tinguished by noting that the oscillating component 
of the light seen by the photomultiplier will suffer 
a phase reversal as the point of observation crosses the 
axis of the plasma column. For odd m, the amplitude 
of the oscillation on axis will be zero. Conversely, if 
a maximum amplitude is found on the axis, the wave is 
either longitudinal, or transverse with even m. These 
conclusions are valid if the perturbations, when repre- 
sented as functions of the radius, show no sign change. 
This is the case for the oscillation discussed in this 
paper. 

A survey of the most important and reproducible odd 
m modes seen in the magnetically confined arc plasma 
in magnetic fields from 0.05 to 0.73 T (1 tesla equals 10 4 
gauss) is presented in figure 1. For convenience in 
our discussion, these modes will be referred to as I, II, 
III, and IV as shown in the figure. Mode II occurs at a 
frequency generally above the ion gyrofrequency, and 
it is this mode that we shall consider in some detail in 
this paper. Modes I, III, and IV appear in an inter- 
mediate field, but I and IV become strongest and 




Figure 1. Transverse, odd m oscillation modes of magnetically 
confined arc plasma. 

Dashed lines: (A) locates the field at which the ion gyroradius is about three times the 
plasma radius; (B) locates the ion gyrofrequency; and (C) locates the field at which the ion 
gyroradius is about one-third the plasma radius. 



sharpest when the conditions for the ideal hydromag- 
netic fluid begin to be satisfied. 

In addition to these odd m modes a persistent even 
m or longitudinal mode is seen at 15 to 20 kHz over a 
wide range of fields. Its frequency and field independ- 
ence suggest that it may be the ion acoustic wave seen 
in a somewhat similar plasma by Alexeff and Neidigh 
[1961, 1963] . There are also others which are neither 
so reproducible nor distinct as the foregoing in the 
region where the ion gyroradius is of the same order of 
magnitude as the plasma column radius. 



TABLE 1. Oscillation and rotation frequencies of plasma column" 
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a All angular velocities in this table (w, W , etc.) are given in megaradians per second. 

h Mode II just begins to appear at this value of ir\L as a band of noise extending from 
0.22 to 0.66 Mrad/sec. oj is taken to be the middle of this band. 

c A second frequency, apparently of another mode, appeared at 0.58 Mrad/sec. It was 
ignored. 



The mode II oscillation behaves like a standing wave, 
in that its frequency is sensitively dependent on the 
length of the discharge between the nozzle and anode. 
(See, for example, columns 1 and 3 of table 1.) This 
was also checked by measuring the longitudinal vari- 
ation of the amplitude of the oscillation. The ampli- 
tude was measured as the rms value, /', of the oscil- 
lating component of the light intensity. /', however, 
contains a dependence on radius and longitudinal posi- 
tion due to inhomogeneity. In order approximately to 
eliminate this R and z dependence, we choose to ex- 
press the amplitude as an equivalent transverse dis- 
placement of the plasma column, AR, defined by: 



AR 



r 



(dl/dR) 



AR is found experimentally to be nearly independent 
of the radius, thus providing evidence that the pertur- 
bation quantities as functions of radius for mode II do 
not show a sign change. The quantity, AR, is shown 
in figure 2A as a function of the z-coordinate. The 
shape of this plot indicates that mode II is a standing 
wave, and that the wavelength is approximately twice 
the column length. 

Some idea of the absolute magnitude of the ampli- 
tude can be obtained by comparing the size of AR, as 
shown in figure 2A with the width of the plasma column 
at half intensity, which is shown in figure 2C in the 
same length units used in figure 2A. AR is thus found 
even at its maximum, near the midpoint of the column, 
to be less than 10 percent of the column width. 
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FIGURE 2. z-dependence of some properties of the magnetically 
confined arc plasma. 

(A). Amplitude of the mode II oscillation, measured in arbitrary units. 

(B). Light intensity on the plasma column axis relative to that at the cathode end. 

(C). Half-intensity full width of the plasma column. Units are the same as in figure 
2A (Three divisions on this arbitrary scale represent approximately 4 mm in half intensity 
full width of the plasma column as seen in the light of the second argon spectrum.) 

(D). Phase angle of the oscillation. This is the same as the angle of the helix described 
by the perturbation for this case, since |m| = 1. 



The departure of the amplitude in figure 2A from 
sinusoidal shape is not at all surprising in view of the 
longitudinal inhomogeneity of the plasma column. 
This is dramatized in figure 2B which shows the lumi- 
nous intensity along the column. The decrease in 
density and temperature which is responsible for this 
is characteristic of d-c discharges of this kind. 

Figure 2D shows that there is a twist in the standing 
wave of the mode II oscillation. If the phase of the 
traveling wave is given by e i{(t>ukz+md \ one would ex- 
pect that of the reflected wave to be e i ^ t ~ kz+md \ 
Composition of such waves gives cos kz e i(<<1>t+m6 \ 
which has a z-dependence of the amplitude similar to 
what is observed, and a perturbation traveling uni- 
formly in the ^-direction, as is observed, but no helicity. 
We have not examined the details of boundary con- 
ditions for reflection of this wave, but have concluded 
that a small attenuation of the wave (consistent with 
the observed width of the resonance) could cause the 
helicity seen in figure 2D. 

That mode II consists of a purely odd m transverse 
wave is seen from figure 3, which shows the amplitude 
of the oscillating component of the light output as the 
photomultiplier traverses the plasma column. The 



Figure 3. Variation of amplitude of mode II oscillation across 
plasma column. 

Reference mark on axis of abscissae locates column axis. 

sharp dip, located at the column axis, precludes the 
possibility of even m transverse or longitudinal mode 
contribution to this oscillation. 

Taking the direction of the magnetic field to be 
positive, the perturbation is found to rotate in the posi- 
tive ^-direction. If the phase factor of the wave is 
e i(a)t+kz+md)^ m mus t De negative. Comparison of 
phases of signals from photomultipliers displaced by 
angles 7r/2 around the column, and of signals from a 
pair of Langmuir probes which subtended an angle of 
20 deg. at the column axis showed that m = — 1 . 

Having thus characterized the mode II oscillation for 
one length (0.54 m) and one field (0.133 T) it remained 
to explore it over as wide a range of physical conditions 
as was practicable. Mode II as seen under the condi- 
tions of figure 1, however, does not represent a good 
starting point. As seen in figure 2B, an attempt to 
vary conditions by changing discharge length may be 
expected to cause the longitudinal inhomogeneity, as 
well as k, to vary. The result is a dependence of fre- 
quency on length that is too complicated to interpret. 
At the same time, an increase in magnetic field (which, 
by decreasing transverse diffusion, improves the longi- 
tudinal homogeneity) is limited, because mode II 
becomes noisy and is lost as its frequency approaches 
the ion gyrofrequency. In order to give maximum 
longitudinal homogeneity, the approach chosen was to 
increase B just beyond the range of mode II, then de- 
crease length until mode II reappeared, and scan it as 
a function of arc length. The data obtained are given 
in table 1, which also includes two typical examples 
of the data obtained at lower fields. 

It was also necessary to determine the angular 
velocity of fluid rotation of the plasma (which occurs 
under the combined influence of the radial components 
of the electric field and the pressure gradient). This 
was measured via the Doppler shift of several spectral 
lines of All (3520.00, 3521.27, 3535.33, 3545.58, 
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3545.84, and 3548.51 A.) Two spectra, one from a 
measured distance either side of the plasma column 
axis, were photographed side by side in the third order 
of a 1.5 m. Czerny-Turner spectrograph with a 
1200 line/mm grating (Applied Research Laboratories 
Quantograph, Model QSMP. The dispersion of this 
instrument is 1.52 A/mm in the third order.) The 
photographic plate was not moved — the spectra were 
positioned by masking the slit. After removing the 
mask, a mercury spectrum was photographed over the 
two argon spectra for reference. Three lines of this 
reference spectrum (2652.04, 2653.68, and 2655.12 A), 
when photographed in fourth order, were conveniently 
located for measuring the shifts between the members 
of each pair of argon spectra. The shifts, which varied 
from 1.7 to 9.3 /x (each value a mean of six determina- 
tions) were measured with a photoelectric scanning 
comparator which could reliably reproduce settings 
on a spectral line to 1 /x. The angular velocities, com- 
puted from the measured Doppler shifts and the radius 
at which the spectra were taken, are recorded in 
column 4 of table 1. 

3. Density and Temperature Distributions 

The density distributions were measured by means 
of the ion saturation current to a Langmuir probe, these 
having been shown to be reliable for this discharge 
in previous work [McLane, Takeda, Thomas, and 
Thompson, 1965]. For convenience in using the densi- 
ties in the calculations to follow, they were graphi- 
cally fit with Gaussian profiles from the axis out to 
5 mm radius. The constants for this fit: 



Tie = Tied exp (— a n R 2 ) 



(1) 



are recorded in columns 11 and 12 of table 1. The 
scatter of the points about the Gaussian was essentially 
random and varied from ± 3.5 percent for the longer 
arc to ± 16 percent for the shortest ones, for which the 
density gradient was larger and the positioning error 
of the probe was more serious. 

The electron temperatures were taken from the 
logarithmic slope of the electron collection portion 
of the Langmuir probe characteristics. These, too, 
were fit by Gaussian profiles, and the constants from 



T e = T eo exp (— a T R 2 ) 



(2) 



are recorded in columns 9 and 10 of table 1. For the 
long arcs, the scatter of points about the Gaussians 
is ± 5 percent if we omit several points near the axis, 
which appear to be too high (possibly due to over- 
heating of the probe in this very flat temperature pro- 
file). This seems justified for the longest arcs, 
(7r/L = 5.8 m~ l ) because the conductance of the plasma 
having the profile obtained from the probe data by 
neglecting these high points agrees to within 10 percent 
(6% for T e ) with that measured with a small radio 
frequency current superposed on the 8 A d-c arc 
operating current. (The method of measurement will 
be reported elsewhere.) For the short arcs, there 
appears to be no problem in the probe measurements 



other than positional inaccuracy and the scatter of 
all points is ± 8.5 percent. This is substantiated in 
one case (tt/L = 15.5 m" 1 - 5 = 0.186 T) by the RF 
measurement, which gave a conductance within 15 
percent (10% for T e ) of that calculated from the probe 
data. 

4. Comparison of Experiment With Theory 

For an inhomogeneous, rotating, cylindrical plasma 
column, Weinstock [1964] has derived the following 
dispersion equation for a generalized hydromagnetic 
wave, including the effect of compressibility and finite 
ion gyrofrequency: 

ay = {I -\m\W-\m\W c l2-k 2 v A 2 1 o) bi 

+ [(W+ \m\WJ2 + k 2 v A 2 lmi) 2 + 2k*v A 2 ]* (3) 

wherein: 

co is the angular frequency of the wave, 

W is the angular frequency of fluid rotation of the 

plasma, 
W c is the angular frequency associated with the dia- 
magnetic current, 
k is the wave number, 
va is the Alfven velocity, and 
o)bi is the angular ion gyrofrequency. 

W, W c and va are to be determined at the axis of the 
column. W c is given by 



Ben e fi RB 



(4) 



wherein p e is the electron pressure, E is the electric 
field, and the other symbols have their previously 
defined, or usual, significance. The Alfven speed is, 
of course, BK/jup)*. In writing (3) and defining 
Vai small correction terms given by Weinstock to 
account for the effect of the arc current are omitted. 
These are negligible for the plasma used in this work. 
It is useful to notice that if the terms W, and W c , 
which are introduced by the inhomogeneity of the 
plasma, are set equal to zero, (3) reduces to: 

CO =-k*VA 2 l«>M + [(k*VA 2 /(O bi ) 2 + 2k*V A ^ 2 ( 5 ) 

Aside from factors of 2, which are introduced into 
(3), and hence, (5), by the finite cylindrical geometry 
of the plasma, (5) is equivalent to the dispersion 
equation for the slow hydromagnetic wave in the 
infinite homogeneous plasma for frequencies much 
lower than the ion plasma frequency. Further, if 
com >> kvA, (5) reduces to 



co = 2 1 l 2 kv A , 



(6) 



which, aside from the 2 1/2 resulting from geometry, 
is the Alfven wave. 

In (4), the second term requires a knowledge of the 
gradient of the plasma potential. Since this is of the 
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order of the ion temperature divided by the plasma 
diameter and the electron temperature is appreciable, 
it is difficult if not impossible to measure. However, 
on the basis of the estimated ion temperature, this 
term should be small — of t he order of 10—15 percent 
of W c under all conditions reported in the table. 
It was therefore omitted in calculating the values given 
in column 7 from the density and temperature profiles. 
It is then straightforward to solve (3) for JV, and 
compute the values given in column 5 of the table 
from the data presented in the other columns. Com- 
parison of the measured and calculated fluid rotation 
velocities shows that: 

(1) The agreement is poor for the data given in the 
first two and the last lines of the table. 

(2) In the other cases, viz, 7r/L values of 6.9 to 15.5 
m _1 in a magnetic field of 0.186 T, the agreement 
appears to be good. 

The effect of the errors estimated in the previous 
section for the density and temperature measurements 
is easily taken into account by solving (3) for W and 
differentiating it with respect to appropriately chosen 
variables. It is found that W is quite insensitive to 
W v — and therefore its value, as calculated from 
a) and other experimental quantities, is insensitive to 
experimental errors in the determination of the density 
and temperature profiles in all cases recorded in table 
1. Further support to this conclusion is provided by 
figure 4, which shows that variation of pressure and 
discharge current over rather wide ranges (which have 
a strong effect on density and temperature distribution) 
have a much smaller effect on the frequency of the 
mode II oscillation than does variation in magnetic 
field (fig. 1). 

The major contributors to the errors in W are the 
wave number, A:, the Alfven speed, va, and especially 
the frequency, co, whose absolute error appears al- 
most directly in the calculated W . If no allowance is 
made for the frequency error, since a) and W (observed) 
were measured simultaneously, allowances are made 
for density and temperature errors of the order esti- 
mated in the previous section, and an error of 10 per- 
cent is assumed in setting k equal to 7r/L, we can come 
to the following conclusions: 

(1) Differences between W obs and W calc of the 
order seen in the first two lines and the last line of 
the table can probably not be accounted for by any 
reasonable estimate of experimental errors. 

(2) Differences of the order seen in all the rest of 
the table (lines 3 through 8) can readily be accounted 
for in terms of experimental errors. 

The data in the first line of the table were taken at 
a comparatively low field and long length. The effect 
of the longitudinal inhomogeneity shown in figure 2B 
is therefore expected to be largest in this case. It 
has already been stated that the ^-dependence of the 
experimental frequencies is complicated here. In 
addition, we cannot be sure that the quantities W c , k, 
and Va, based on measurements made in the midplane 
of the arc column, are the proper mean values to use 
in order to take account of the longitudinal inhomo- 
geneity. The calculated values of W may therefore 
be incorrect. 
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Figure 4. (A). Dependence of mode II frequency on pressure 
(Discharge current, 8A; magnetic field, 0.106 T; length, 0.54 m). 
(B). Dependence of mode II frequency on discharge current 
(pressure, 0.95 m torr.; magnetic field, 0.106 T; length, 0.54 m). 

The results for the largest values of tt/L (lines 2 
and 9) are best viewed in terms of the limiting form of 
(3) for large A^vJ/cww 



o,' S a»-(2-W)r*a,„-|^ 



(7) 



For \m\ — 1, this sets an upper limit on the frequency 
measured relative to the rotating plasma which is 
equal to the ion gyrofrequency. In the sequence of 
lines 3 to 9 of the table, it is clear that when co' clearly 
surpasses co/>; with increasing A 2 t^/o>w, the measured 
and calculated values of W no longer agree. This 
is not surprising, for this is just the condition under 
which the wave described by (3) can no longer exist — 
transition to another k or m value will not, in general, 
lead to a possible solution if a/ > o)bu While the 
plasma nevertheless appears to continue oscillating 
with a similar frequency in a transverse mode above 
this resonance, the detailed study of the structure of 
the wave which would be necessary to identify the 
mode is difficult for such short discharges, and has 
not been done. 

5. Conclusion 

A transverse, standing wave oscillation has been 
seen on the plasma column of a magnetically confined 
arc. Over a two to one range of wave numbers and 
frequencies, it obeys a dispersion equation applicable 
to a generalized hydromagnetic wave in a rotating, 
inhomogeneous, cylindrical plasma. The agreement 
is lost when the frequency relative to the rotating 
plasma surpasses co&i, and similar waves in which the 
fluid rotation velocity is abnormally small occur. 
These have not been identified. 
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